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Assuming that curvature perturbations and gravitational waves originally arise from vacuum 
fluctuations in a matter-dominated phase of contraction, we study the dynamics of the cosmological 
perturbations evolving through a nonsingular bouncing phase described by a generic single scalar 
held Lagrangian minimally coupled to Einstein gravity. In order for such a model to be consistent 
with the current upper limits on the tensor-to-scalar ratio, there must be an enhancement of the 
curvature huctuations during the bounce phase. We show that, while it remains possible to enlarge 
the amplitude of curvature perturbations due to the nontrivial background evolution, this growth 
is very limited because of the conservation of curvature perturbations on super-Hubble scales. We 
further perform a general analysis of the evolution of primordial non-Gaussianities through the 
bounce phase. By studying the general form of the bispectrum we show that the non-Gaussianity 
parameter /nl (which is of order unity before the bounce phase) is enhanced during the bounce phase 
if the curvature fluctuations grow. Hence, in such nonsingular bounce models with matter given by a 
single scalar field, there appears to be a tension between obtaining a small enough tensor-to-scalar 
ratio and not obtaining a value of /nl in excess of the current upper bounds. This conclusion may 
be considered as a “no-go” theorem that rules out any single field matter bounce cosmology starting 
with vacuum initial conditions for the fluctuations. 
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I. INTRODUCTION 

As was realized in [Ilia, there is a duality between the evolution of curvature fluctuations in an exponentially 
expanding universe and in a contracting universe with the equation of state of matter. In both cases, curvature 
fluctuations which originate as quantum vacuum perturbations on sub-Hubble scales acquire a scale-invariant spectrum 
at later times on super-Hubble scales. The observed small red tilt of the spectrum of curvature perturbations which 
has now been confirmed by observations (see e.g. mm can be obtained in an expanding universe by a slow decrease 
of the Hubble constant during the period of inflation [^, whereas in a matter-dominated phase of contraction a small 
cosmological constant (with magnitude comparable to what is needed to explain today’s dark energy) yields the same 
tilt (see alternatively 0 )- To avoid reaching a singularity at the end of the contracting phase, it is necessary to 
either modify gravity or consider matter violating the null energy condition (NEC). Then it is possible to obtain 
nonsingular bouncing cosmologies which have the potential to yield an explanation for the structures in the universe 
which we now observe. This scenario of structure formation alternative to inflation is called the “matter bounce” 
scenario (see e.g. EllS] for reviews). 

Examples of modified gravity models which yield bouncing cosmologies include the “nonsingular Universe” construc¬ 
tion of nniiii], nonlocal gravity actions like |12j , or Hofava-Lifshitz gravity |I3j . It is in general very hard to study 
the evolution of fluctuations in these models. We will hence focus on models in which the bounce is obtained from 
the matter sector. One method of obtaining a nonsingular bounce with a single scalar field involves the formation of 
a ghost condensate during the bounce phase (see [MUI^ for initial developments). A general problem for bouncing 
cosmologies is the Belinsky-Khalatnikov-Lifshitz (BKL) instability |5D], the fact that the energy density in the form 
of anisotropies will explode and destroy the homogeneous bounce m- This problem can be “solved” by endowing 
the scalar matter field with a negative potential which leads to an Ekpyrotic phase of contraction before the bounce 
[221 - 124] and hence can mitigate the anisotropy problem [25] 

In the matter bounce scenario, primordial quantum fluctuations exit the Hubble horizon while the universe is in a 
matter-dominated contracting phase and the resulting power spectrum of curvature perturbations is scale-invariant 
[DEj. On the other hand, the gravitational wave mode obeys the same equation of motion on super-Hubble scales as 
the curvature perturbations (considering the canonical variables in each case). Hence, before the bounce phase the 
tensor-to-scalar ratio r would be of order unity. Thus, if the perturbations passed through the nonsingular bounce 
unchanged, it would imply that curvature perturbations and primordial gravitational waves would have the same 
amplitude after the bounce. In terms of the tensor-to-scalar ratio, it would mean that r ~ C(l), well above current 
observational bounds I1S1E3. 


^ Such a negative potential may arise from the standard model Higgs field since, based on the recent Higgs and top quark mass measurements, 
the standard model Higgs develops an instability at large field values (in the absence of new physics) ESj. 
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Since the curvature fluctuations couple nontrivially to matter during the bounce phase, whereas the tensor perturba¬ 
tions are determined simply by the evolution of the scale factor a{t), one may expect that the curvature perturbations 
would be enhanced relative to the tensor modes during the bounce. In fact, early calculations indicated that curvature 
perturbations grew exponentially during the bounce phase, hence suppressing the tensor-to-scalar ratio [531 [55] ■ A 
more recent study |29j numerically explored the evolution of scalar fluctuations through a nonsingular bounce model 
similar to the one studied in [22] and found no enhancement of curvature perturbations through the bounce. In light 
of the relevance of a possible enhancement of the curvature fluctuations for the predicted value of the tensor-to-scalar 
ratio, the growth of curvature fluctuations during a nonsingular bounce needs to be reconsidered. This is what we aim 
to do in this paper. 

The second goal of this paper is to carefully track the evolution of the three-point function (bispectrum) of curvature 
perturbations through the bounce. In earlier work |30j it was shown that the bispectrum of curvature fluctuations 
before the bounce phase has an amplitude of the order /nl 0{1) with a specific shape. As we argued above, if 
the perturbations were to pass through the nonsingular bounce unchanged, it would imply a large tensor-to-scalar 
ratio in excess of the observational bounds. On the other hand, if curvature perturbations were to experience a 
nontrivial growth through the bounce, one should expect additional nonzero contributions to the bispectrum coming 
from the bounce phase, and there would then be the danger that the final amplitude of the bispectrum exceeds the 
observational upper bounds from [311 132 j . Thus, a potential conflict looms: either the tensor-to-scalar ratio is too 
large, or else the non-Gaussianities exceed observational bounds. This problem has indeed already been found in a 
model of a nonsingular bouncing cosmology in which a nonvanishing positive spatial curvature is responsible for the 
bounce [221 [21] ■ We will study this issue in the context of the more realistic models in which the nonsingular bounce is 
generated by the matter sector. In particular, we will explore the question in the context of a ghost-condensate bounce. 

We will indeed demonstrate that, at least in our model, the evolution of the curvature perturbations in the bounce 
phase connects the value of the tensor-to-scalar ratio with the amplitude of non-Gaussianities. The suppression of the 
tensor-to-scalar ratio to restore compatibility with the observational bounds requires an enhancement of the curvature 
fluctuations during the bounce phase. Such an enhancement will increase the magnitude of the non-Gaussianities 
to a level inconsistent with the observational bounds on the amplitude of the bispectrum. Based on our result we 
conjecture that there exists a “no-go” theorem in single field nonsingular matter bounce cosmologies which relates the 
tensor-to-scalar ratio and non-Gaussianities, preventing these models to satisfy the current observational bounds. A 
tensor-to-scalar ratio below current observational bounds would imply a too large amplitude of non-Gaussianities, 
whereas non-Gaussianities of order /nl 0{\) would imply a too large amplitude of the primordial gravitational wave 
spectrum. Therefore, a single field nonsingular matter bounce cannot be made consistent with current observations if 
the primordial perturbations arise from vacuum initial conditions. 

Our analysis assumes that both curvature perturbations and gravitational waves originate as quantum vacuum 
fluctuations in the initial phase of contraction. A model with thermal fluctuations (as obtained for example in the 
context of string gas cosmology [23125]) will easily avoid our “no-go” theorem. As shown in [HTIHO] . we obtain 
a tensor-to-scalar ratio much smaller than order unity while obtaining non-Gaussianities which are negligible on 
cosmological scales [41] . 

The paper is organized as follows. We first start with a short review of cosmological perturbation theory in Sec. 
[Ill We then motivate the idea of the no-go theorem proposed in this paper in Sec. |III| In Sec. IV we briefly review 
the general picture of bouncing cosmology in terms of a single scalar field of Galileon type. After that, in Sec. [^ 
we analyze the perturbation equation for primordial curvature perturbations at linear order during the nonsingular 
bouncing phase. We point out under which conditions there can be an enhancement of their amplitude. Then in 
Sec. |V^ we perform a detailed analysis of the bispectrum generated in the bouncing phase of our specific model. We 
combine the analyses of scalar and tensor perturbations together with non-Gaussianities in Sec. |VII1 and we show how 
current observational bounds severely constrain the parameter space of the single field bouncing model. The analysis 
is expected to hold quite generally for single field matter bounce cosmologies. We conclude with a discussion in Sec. 
VIII[ Throughout this paper, we adopt the mostly minus convention for the metric and define the reduced Planck 


mass as M/ = I/SttGn where Gn is Newton’s gravitational constant. 


II. A BRIEF REVIEW OF COSMOLOGICAL PERTURBATION THEORY 

Linear perturbations of the metric about a homogeneous and isotropic background space-time can be decomposed 
into scalar, vector, and tensor modes (see |42j for a review of the theory of cosmological perturbations and |43j for an 
introductory overview). The scalar modes are those which couple to matter energy density and pressure perturbations. 
We call these the cosmological perturbations. Tensor modes exist in the absence of matter - they correspond to 
gravitational waves. In the case of matter without anisotropic stress at linear order in the amplitude of the fluctuations, 
there is only one physical degree of freedom for the scalar fluctuations. For the purpose of computations it is often 
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convenient to work in the conformal Newtonian gauge (coordinate system) in which the perturbed metric for scalar 
modes reads 


= 0 ^( 77 ) ([1 + 2 <i>(ry, a;)] — [1 — 2 <i)(ry, a:)] , ( 1 ) 

where rj denotes conformal time, 0(77) is the cosmological scale factor, x represents comoving spatial coordinates, and 
$ denotes the gravitational potential. For tensor modes, the perturbed metric reads 

ds^ = a^{r]) {drf — [Sij + hij{r], £)] da;*dx^) , ( 2 ) 

where hij is trace-free and divergenceless. 

Let us consider the matter content to be described by a single scalar field of canonical form with Lagrangian density 

Cm = - V{^) . ( 3 ) 

Note that we take the scalar field to be dimensionless throughout this paper as a convention. Linear perturbations of 
the scalar field then have the form 


Hv, x) = (j)o iv) + Hiv, x), ( 4 ) 

where (j)o is the unperturbed homogeneous part of 4>. In the scalar sector, metric and matter perturbations couple to 
one another, so it is useful to define a linear combination of these perturbations, 

7 ^ = + $ . ( 5 ) 

Vo 

There are two reasons for focusing on this variable. First of all, it gives the curvature fluctuation in comoving 
coordinates (coordinates in which the matter field is uniform), and is hence the variable we are interested in computing. 
Second, it is simply related to the Sasaki-Mukhanov [331 |3S] variable v in terms of which the action for cosmological 
perturbations has canonical form. Note that in the above, T-L = a! ja is the conformal Hubble parameter and a prime 
denotes a derivative with respect to conformal time. In fact, the Sasaki-Mukhanov variable is 


V = zK , ( 6 ) 

with 

z = a^Mp. ( 7 ) 

The equation of motion that results from expanding the perturbed action for gravity and matter to second order is 
given by 

v'k + “ t ) ^ ^ 

The equation is written in Fourier space, where k represents the comoving wave number of the curvature perturbations, 
and Cs is the speed of sound which is equal to one for a scalar field with canonical action ([^. Similarly, for tensor 
modes the Mukhanov variable is 


^ = ah, 


( 9 ) 


where h is the amplitude of the polarization tensor hij (the two polarization states evolve independently at linear 
order and obey the same equation of motion) and the resulting equation of motion is 


hk 


Csk^ -Tifc = 0 . 

a 


( 10 ) 


Alternatively, without the use of the Mukhanov variables, the equation of motion for curvature and tensor perturbations 
can be written as 


z' 

7^fc -b 2-7^; -b c^.k^Tlk = 0 , (11) 

Z 

h'k + 2 h'j^ -b c^k'^hk = 0 , 
a 


(12) 
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respectively. 

Finally, let us introduce the scalar perturbation variable 

d = $ -I- 

^ “ sn{i + w) ’ 


(13) 


where w = P/p is the equation of state parameter (P is the pressure and p is the energy density). On super-Hubble 
scales, i.e. for fc <C P, this variable is equivalent to the curvature perturbation variable TZk |46j . In other words, 
Ti-k = Cfci thus, throughout the rest of this paper, we will use TZu and C,k interchangeably to denote curvature 
perturbations on super-Hubble scales. 


III. OUTLINE OF THE NO-GO CONJECTURE 

As explained in the introduction, a careful study of the evolution of curvature perturbations and the production 
of non-Gaussianities during a nonsingular bounce may lead to a “no-go” theorem, the impossibility of obtaining a 
sufficiently small tensor-to-scalar ratio while maintaining a bispectrum with an amplitude smaller than the current 
observational bounds. In this section we will provide a qualitative analysis of this problem by giving simple estimates 
of the tensor-to-scalar ratio and of the amplitude of the bispectrum assuming that the curvature fluctuations undergo 
some growth through the bounce phase. We first start by setting up the matter bounce formalism. 


A. Fluctuations in the matter bounce 

In the matter bounce, primordial quantum fluctuations originate on sub-Hubble scales during a matter-dominated 
contracting phase and exit the Hubble radius during this phase. The perturbations then remain on super-Hubble 
scales as the universe contracts and passes through the bounce phase, except for a very small time interval right at the 
bounce point (at which time the Hubble radius goes to infinity). The ffuctuations with wavelength of cosmological 
interest today will then reenter the Hubble radius in the standard radiation or matter-dominated expanding phases. If 
the bounce is completely symmetric, then fluctuations which exit the Hubble radius in the matter phase of contraction 
reenter the Hubble radius in the matter phase of expansion. However, we expect the bounce to be asymmetric and 
entropy to be generated during the bounce. In this case, the radiation phase of expansion is longer than the radiation 
phase of contraction. 

To understand the evolution of quantum fluctuations in a contracting universe, one needs to determine the form of 
the variable z and then solve Eq. ([^. Using the Friedmann equations, the time derivative of the Hubble parameter is 
given by 


H = - 



(14) 


where a dot denotes a derivative with respect to cosmic time, t, and the subscript 0 indicates that we are referring to 
the background field. Defining the parameter e, 


€ = — 


H 

IP ’ 


(15) 


and using Eq. (14), one finds 


z = a^Mp = aV^Mp . 
H 


(16) 


It is straightforward to show from the Friedmann equations that 


e = -(1 + w) 


(17) 


so for a matter-dominated contracting universe with ru = 0, we have e = 3/2. As a consequence, z = a\/3Mp and 


z 

z 


a 

a 


(18) 
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and we conclude that the scalar and tensor fluctuations evolve in exactly the same way. This is not true in general 
since w can vary in time. For example, in the case of inflationary cosmology, we recognize e as the slow-roll parameter 
and it is time-dependent. 

In a matter-dominated contracting universe, the scale factor scales as a ~ (—~ 77 ^, and since cf = 1 for a 
canonical scalar field, the equation for the Sasaki-Mukhanov variable is 

v'k + = 0 ■ ( 19 ) 

On super-Hubble scales, the term is negligible, and so the solution reads 

Vkiv) = CiTf+C 2 T]~^ ■ ( 20 ) 

Using the fact that Vk = z(^k, the first term yields C^k ~ constant, but in a contracting universe, the second term is the 
dominant solution, 

Cfc ~ , (21) 

which implies that curvature perturbations grow in a contracting universe. In fact, the growth rate is precisely the 
correct one to convert an initial vacuum spectrum into a scale-invariant one (see e.g. [9] for a review). 


B. Bound from the tensor-to-scalar ratio 


The tensor-to-scalar ratio is defined as 

Vt{K) 


( 22 ) 


where fc* is the pivot scale which is used to parametrize the power spectra for tensor and curvature perturbations. The 
individual power spectra are defined by [37] 


Vt{k) = 2Vh{k) = 2 X 16n^\hk\^ = 

ZTT^ TT^ 

V (k) = —\r P = 


(23) 

(24) 


respectively. The factor of 2 in the first step of the first line comes from the two polarization states of gravitons and 
the factor of IGtt is a convention reflecting the fact that it is IGirMph which yields the canonical action of a free scalar 
field in an expanding background [37] . 

As we found in the previous subsection, 2 = a'/SMp for the matter bounce, so the scalar power spectrum becomes 


V^k) 


k^ kfcP 

67r2 a2M2 ’ 


and furthermore, the tensor-to-scalar ratio becomes 


r = QGtt 


Mfe. 


Vk, 


2 


Mi 


(25) 


(26) 


where the factor reflects the fact that we have defined Vk to have dimensions of mass, whereas ^k is dimensionless. 

Since z" jz = o!' ja for the matter bounce, the evolution of scalar and tensor modes given by Eqs. Q and (10), 
respectively, will be identical. In addition, if they originate from the same quantum vacuum, then Vk{i]) = Mp^kiv)- 
Consequently, we find that r = 967r. If perturbations passed through the bounce unchanged, it would result in r = 967r 
at the beginning of the standard big bang cosmology phase which is three orders of magnitude larger than the current 
observational upper bound. 

To gain some intuition on the effect of passing through the bounce phase, let us assume that curvature perturbations 
are enhanced by an amount A(^fc through the bounce, i.e. 


(k{r]B+) = CkiVB-) + 


(27) 
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where f}B± denote the conformal time before (—) and after (+) the bounce. Then, the tensor-to-scalar ratio measured 
after the bounce becomes 


r['qB+) = 967r 


hkAVB+) 

CkAVE-) + AC/c. 


(28) 


Assuming that tensor modes remain constant through the bounce, i.e. hkAjB-) = hk(j]B+), one finds that 


1 + 


CkAVB-) 


Avb-) 

'r{'nB+) 


(29) 


Taking the value of the tensor-to-scalar ratio before the bounce to be what we found earlier, i.e. r{r]B-) = 967r, and 
demanding that the tensor-to-scalar ratio is sufficiently suppressed after the bounce so that it satisfies the observational 
bound r{r]B+) < 0.12 (95% CL from [23|48j), we find that curvature perturbations must be sufficiently enhanced 
during the bounce phase so that 


or using the triangle inequality. 


1-f 


Aq. 

CkAvB-) 


A 50.1 , 


Aq. 

CkAvE-) 


A 49.1 . 


(30) 


(31) 


C. Bound from the bispectrum 


The primordial bispectrum, B(^, is defined in terms of the three-point function as 

(C(fci)C(fc2)C(4)) = {27rAS^^Aki + k2 + h)BAki,k2,k3) , 

which we can rewrite as 

(C(^i)C(fc2)C(4)) = (27r)75(3) 






n, 


j^A{ki,k2,k3) 


(32) 


(33) 


where ki = \ki\ and where the index i runs from 1 to 3. The function A{ki, ^ 2 , fcs) is known as the shape function and 
its amplitude defines the nonlinear parameter /nl via 


hL{ki,k2,k3) 


10 A{ki, k2, ks) 

^ E.kf 


(34) 


Of particular interest is the local form of non-Gaussianities for which one of the three modes exits the Hubble radius 
much earlier than the other two, i.e. ki k 2 = k^. For this case, one can write 

C{x) = Cgix) + ^/Nr^Cg(^)^ , (35) 

where q is the Gaussian part of q 

In order to compute /nL) one must evaluate the three-point function. To leading order in the interaction coupling 
constant, the three-point function is related to the interaction Lagrangian, Lint, via |49j 


(C(i,fci)C(Cfc2)C(Cfc3)) =ij^ di ([C(t,fci)C(i,fe)C(^, 4), Aint(t)]) , (36) 

where the square brackets denote the commutator and where U denotes the initial time before which there is no 
non-Gaussianity. The interaction Lagrangian is obtained by evaluating the action up to third order in perturbation 
theory 


Aint(t) = J d^x C3{t,x) , 


(37) 

















(38) 


and for a canonical scalar field, the Lagrangian density for C to cubic order is given by [49] 


M2 

/(C) 


if ~ i) + «'<(^C)" - 2e^a^ad0{dx) + + /(C)^ , 

- w<C-^d.(d.x + ^d-%dA9.xdK), 


(39) 


where i9 ^ is the inverse Laplacian and where we define x = 9 ^C- Also, the equation of motion for ( coming from the 
second order perturbed Lagrangian density £2 is given by 


S£2 

<5C 


C2z2 


4(a2:^C) - 

at a 


d\. 


(40) 


As we saw in Sec. |III A] curvature perturbations grow on super-Hubble scales during the matter-dominated contracting 
phase until the bounce phase. While on super-Hubble scales the spatial gradient terms are negligible, i.e. di(, diX — 0, 
the growth in ^ implies that the interaction Lagrangian is dominated by 


M2 



1 

H 


c4(-^C) . 


(41) 


As was first shown in |30j . the production of non-Gaussianities on a comoving scale k is dominated by the period 
between when the scale crosses the Hubble radius in the phase of matter contraction until the onset of the bounce 
phase, and the resulting non-Gaussianities are of order /nl ~ C>(1). For example, for the local shape, the authors of 
[30] found = -35/16. 

Following what was done in the previous subsection, let us now assume that curvature perturbations grow during 
the bounce phase. For simplicity, let us assume that they grow linearly in time with constant rate 


c = 


Ats ’ 


(42) 


where the duration of the bounce is given by Ats = tB+ — ts-- Then, in the limit A: —>■ 0 on super-Hubble scales, the 
contribution to the three-point function coming from the bounce phase is schematically given by 


{atB+f) 


bounce 


C(^B+)^ 

Mp 





dt a{t)^ 


' ts- 



2 


((ts-) + 


Ats 


(t-tB-) , 


(43) 


and one expects that the dominant contribution to /nl that results from evaluating the three-point function would 
scale as 


/nl 


AtB ^ 


(44) 


plus terms of order which would be subdominant for a large amplification A/. 

We already see that a growth in the curvature perturbations during the bounce, A((, would enhance /nl- From 
the previous subsection, we expect A/' to have a lower bound to match current observational bounds on r, and thus, 
we expect to find a lower bound on the amount of non-Gaussianities that are produced during the bounce phase. 
However, we cannot determine whether this contribution will be significant to /nl C>(1) and whether the resulting 
lower bound will exceed current observational bounds without going into the details of the calculation. 


D. The no-go theorem 


Now, let us state our conjecture. 

Conjecture 1 For quantum fluctuations originating from a matter-dominated contracting universe, an upper bound 
on the tensor-to-scalar ratio (r) is equivalent to a lower bound on the amplification of curvature perturbations (A//// 
which in turn is equivalent to a lower bound on the amount of primordial non-Gaussianities (/nl/- Furthermore, if 
the initial quantum vacuum is a canonical Bunch-Davies vacuum with Cg = 1, if the nonsingular bounce phase is due 
to a single NEC violating scalar field, and if general relativity holds at all energy scales, then satisfying the current 
observational upper bound on the tensor-to-scalar ratio cannot be done without contradicting the current observational 
upper hounds on /nl (and vice-versa). 

In the rest of this paper, we will give an example of realization of this conjecture. 
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IV. A BRIEF REVIEW OF SINGLE FIELD BOUNCING COSMOLOGY 


In the context of Einstein gravity, matter which violates the null energy condition must be introduced in order to 
obtain a cosmological bounce. A simple toy model is quintom cosmology, i.e. a model in which a scalar field with 
opposite sign in the action compared to a usual scalar field is introduced, and it is arranged that this field comes to 
dominate late in the contracting phase, thus yielding a nonsingular bounce [50]. A specific realization of this can 
be obtained in the Lee-Wick theory m- These models, however, suffer from a ghost instability |52|- To avoid this 
instability (at least at the perturbative level) one can make use of the ghost condensation mechanism [TO] or the 
Galileon construction [nms] • ^ These mechanisms involve a modified kinetic term in the action. 

As mentioned in the introduction, bouncing models typically also suffer from the anisotropy problem, and to mitigate 
this problem, one can build into the scenario an Ekpyrotic phase of contraction which occurs at some point after 
the matter phase of contraction. Specifically, one can use a single scalar field with a kinetic term designed to yield a 
nonsingular bounce, and a potential energy function with a negative potential over some range of field values which is 
designed to yield Ekpyrotic contraction [22] . In this approach, a second scalar field with canonical kinetic term and 
with quadratic potential can be used to represent the regular matter of the Universe [23] . In this paper we will not 
consider the role which this second scalar field may play (for some ideas see |5B|) but only consider the field (j) which 
generates the Ekpyrotic contraction and the nonsingular bounce. 

Throughout this paper, we assume only Einstein gravity plus matter. Thus, the action is given by 


S = 



^ m2 




(45) 


where g is the determinant of the metric, i? is the Ricci scalar, and Cm is the matter Lagrangian. We assume that the 
matter content is dominated by only one scalar field (<()) before reheating (the energy density of matter created via 
reheating becomes only important after the bounce phase - see [24] 1. Thus, for the dynamics of the matter-dominated 
contracting era and the bounce phase to be described by second order equations of motion, we consider a Lagrangian 
of the most general form |5U] 


Cm — X) -|- X)\I\(j) + C 4 + £5 , 


(46) 


where the kinetic variable X is defined as 

A = , (47) 

and where the d’Alembertian operator is defined as 

. (48) 


We do not write down the explicit form that £4 and £5 can take here, but the key point is that they involve higher 
order derivatives. If we assume that the energy scale at which the bounce occurs is low enough so that higher order 
derivative terms in the Lagrangian are negligible, then we can assume that £ 4 , £5 « 0. 

Eor the bounce to be nonsingular, the above Lagrangian must violate the null energy condition (NEC) at high 
energies. To do so, we assume the first term of the Lagrangian to have the form 


A(</., A) = M2[1 - 5 (<(,)]A + /?A2 - I/(^) 


(49) 


where j3 is some positive constant. We see from Eq. (49) that when g{(j)) > 1, the sign of the kinetic term is reversed 
and a ghost condensate which violates the NEC is formed [HHISKISIIIS]. Eor this reason, one typically chooses the 
function g{<j)) to have the form 




_ 2^0 _ 

_|_ gbgy/2jp<l, 


(50) 


where p and bg are positive constants. As ()>—?. 0 at the bounce point, g{4>) —>■ go, and the constant go is naturally 
chosen to be go > 1 to allow the NEC violation. We can also see from the form of g{(f>) above that as </> goes away from 


2 


Alternative possibilities of alleviating this instability may be achieved by considering various modified gravity implementations such as 
models of extended F(R) gravity |53II54| . modified Gauss-Bonnet gravity |55| . and torsion gravity scenarios |56lI57| . 
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0 and as the kinetic variable X becomes small outside the bounce phase, g{(f>) rapidly goes to 0 and the Lagrangian 
recovers its canonical form. 

The potential V (</>) can be chosen in order to obtain an Ekpyrotic phase of contraction. This can be done by means 
of a potential which is negative for small values of \4>\, but which approaches V = 0 exponentially at large positive and 
negative field values. Specifically, we have chosen the potential 


E(<^) 


2Vo 


(51) 


where Vq, q, and by are positive constants. Without the second term in the denominator, one obtains the potential 
postulated in the Ekpyrotic scenario EO]. 

One can then parametrize the background evolution during the bounce phase as follows. The Hubble parameter 
grows linearly in time, passing through zero at the time t (the bounce point). 



1 

II 

(52) 

where T is a positive constant. 

The scale factor immediately follows. 



a{t) = 

(53) 

Also, the scalar field evolves as 


(54) 


Since qb and ts can be arbitrarily redefined, we see that the parameters which describe the bounce phase are T, rj)B, 
tB- (or tB+ assuming a symmetric bounce), and T. First, T gives the growth rate of the Hubble parameter. Second, 
(j)B gives the maximal growth rate of the scalar field. Third, AIb/T gives the dimensionless duration of the bounce. 
They can be related to the Lagrangian parameters via (see [22 (22] ) 


. /2(go - 1) ,, 

I _ ^ 


(55) 

(56) 


where Hb+ = T{tB+ - tB)- 

Given the model we have discussed and the bounce solution which we have given in parametric form, we will now 
follow the evolution of the curvature fluctuation variable C through the nonsingular bounce phase. 


V. EVOLUTION OF CURVATURE PERTURBATIONS DURING THE BOUNCE 


As we saw in Sec. |Hj the equation of motion for curvature perturbations [Eq. 0] can be written as 


7^ 


TZ'i. + c^„k'^TLk = 0 . 


fc W ^2 k 


For a noncanonical Lagrangian of the form of Eq. (46), the variable 0 and the sound speed are given by 

2M^a^,y^r 


z^ = 


{2M^H - G.x<^ 3)2 ’ 


1 r 


c, = 


V L 


K^x + ^H(t)G^x - 


2 Ml 


- 2G,^ + G,x4><y^ + {2G,x + 


where a comma denotes a partial derivative and where we defined 

3 


V = K^x + (t> K,xx + 


i4gi2 


2 M 2 


G^x + QH4 >G^x + 3Hcj)'^G.xx - 2G,^ - ^G^^x ■ 


(57) 

(58) 

(59) 


(60) 
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FIG. 1: Sketch of the different regimes in the bounce phase (not to scale). The horizontal axis represents physical time. The 
green solid curve shows 2MpH{t) and the dashed version depicts its absolute value. The bell-shaped blue curve represent 
\G.x it)\ 4>^{t), where we take G,v = 7 to be a positive constant for simplicity. Regimes I, II, and III, defined by Eqs. (63l, (64l, 
and (|65[), are depicted by the pink, purple, and cyan regions, respectively. 


As explained in Sec. |III A[ the perturbation modes that are of cosmological interest today were on super-Hubble scales 
during the bounce phase (except in the immediate vicinity of the bounce point), and thus we are most interested in 
the infrared (IR) regime of Eq. ( |57| ). In the limit fc <C 'H, and recalling that TZk and Ck are equivalent quantities in this 
limit, the equation that we want to solve is 


dri 



C' = o, 


(61) 


where we drop the k index when it is clear that we are on super-Hubble scales. It is obvious from the above equation 
that one solution is the constant mode solution, = 0, that one expects on super-Hubble scales, e.g. in inflation 
(see, however, [63]). More generally, the solution to Eq. (61) can be written as 


Civ) = C{v^) 


C^ivi) 

z^iv) 


(62) 


where rji denotes the initial time where the initial conditions are set. The evolution of C is thus governed by the 
evolution of , and we notice from the denominator of Eq. (58) that the evolution of z'^ has different regimes of 
interest: 


Regime I : 

2 Ml\H{t)\ » \G,xit)\4C{t) , 

(63) 

Regime H : 

2 Ml\H{t )\« \G,x{t)\4C{t) , 

(64) 

Regime HI 

: 2 MlH{t) « G,xmHt) ■ 

(65) 


We represent these different regimes in Fig. and we explore the consequences of each regime in the following 
subsections. 


A. Evolution in Regime I 


When Eq. (63) is valid, the expression for z'^ reduces to 


/ !-£/((/)) 

~ m \ 


2 


(66) 
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Since the bounce phase is defined by g{(j)) > 1 and since must be positive to avoid ghost instabilities, the model 
parameters must be chosen such that this regime does not occur during the bounce phase. Outside the bouncing phase, 
the equation of motion in Regime 1 reduces to the standard one. 


B. Evolution in Regime II 


As the bounce point approaches, H{t) goes to zero and we can expect Eq. (64) to be valid. To explore this regime, 
let us simplify the treatment by setting 

G{^,X)=jX (67) 

for some positive constant 7 , so the regime becomes 

2M^\H{t)\ « • ( 68 ) 

Using the parametrizations introduced in the previous section, this condition can be rewritten as 

^ , (69) 

where we defined At = t — ts- Since Atg/T determines the dimensionless duration of the bounce, remaining close to 
the bounce is equivalent to demanding that |At|/T ^ 1. In particular, if we demand that 


|At| < 


^ 7^ 1 

v^’2M2Tj ’ 


(70) 


then it is ensured that we are in the regime set by Eq. ( p^ . Thus, the expression for given in Eq. (58) reduces to 


7 


in this regime. In fact, there exists a time interval, which we define as [tamp-, iamp+] with tamp± = ts ± Atamp, where 
the above approximation for z^(t) is certainly valid. We note that this expression is everywhere finite in that interval, 
so the solution to Eq. (61) can be directly written as 


at) = m 


a{ti)z‘^{ti) 

a{t)z'^{t) 


(72) 


where the initial condition must be taken in the interval, i.e. ti S [tamp-, tamp+]i so logically we take ti = tamp-- Also, 
the solution will only be valid up to tamp+- Inserting Eq. ( |71[ ) and using the parametrizations introduced in the 
previous section, one finds 


C(^) — C(^amp—) T C(tamp—) 
~ C(^amp—) “t” C(tamp—) 
~ C(^amp—) “t” C{^amp—) 


dt 


^(tamp— ) 
a{i) 


0(t) 


0(tamp— ) 


(^(tamp— ) 

as 

^(tamp— ) 

as 


4>b 


erf 


t-ts 

T 


3T2T \ 

2+ —l-erf 


^ 0(tamp— ) J 

^ ^(tamp— ) J 

tamp— tfj 

• -J, ' 


dt exp 


8 + 6T2T 




2 + 


3r2T \ 


(73) 


Close to the bounce point, the scale factor remains nearly constant, so a(t) ~ as- This implies that T(At)^ ^ 2, or in 
other words, that H{t)At ^ 0{1). We will assume this to be valid throughout the rest of this paper whenever we are 
in the time interval |At| < Atamp- Therefore, the solution for ({t) reduces to 


C(^) — C(^amp—) + C(tamp—) 


tV^ 


^(tamp— ) 


erf 


^-^ 02 ] - erf { 


(74) 
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FIG. 2: Sketch of the evolution of curvature perturbation on super-Hubble scales as a function of physical time t (not to 
scale). The beginning of the bounce phase, the bounce point, and the end of the bounce phase are denoted by ts-, ts, and ts+, 
respectively. We defined tamp± = ts ± Atamp, and ts is the time at which ^ oo. The purple region corresponds to Regime II 
of Fig. where C, grows at most linearly. The cyan region corresponds to Regime III of Fig. where C, is almost constant. 


From the above solution, we see the constant mode and the growing mode. Whether the constant or the growing 
mode is dominant depends on many factors. For instance, the duration of this regime and the growth rate will play a 
crucial role. From the properties of the error function, we note that the growing mode grows at most linearly in time. 
Furthermore, the growth rate is maximal at the bounce point ts and it is given by Cmax — C(^amp-)[^B/<(>(^amp-)]^- 


C. Evolution in Regime III 


One can notice from Eq. (62) that if —>■ oo, then C,' —t 0, and curvature perturbations remain constant on 
super-Hubble scales. One can see from Eq. (58) that this happens at some physical time tg (or rjg in conformal time) 
when 


2M^H{ts) = G,xits)<^Htg) . (75) 

At this point, the equation of motion for the curvature perturbations becomes singular, and furthermore, the Mukhanov 
variable Vk = zTZk diverges. For this reason, the evolution of the curvature perturbations has been explored in another 
gauge, the harmonic gauge (first introduced in the context of cosmological perturbation theory in |64j). where this 
singularity may disappear. Using the harmonic gauge, it has been shown in [29] that at 77 ^, 


dTZk 

drj 


= 0 


(76) 


for all k modes. Carefully dealing with the singular equation of motion in the conformal Newtonian gauge, one can 
find that in the IR limit, the solution in conformal time close to the singular time r]g is (see Appendix [A]) 


Civ) = Civ^) + C'iv^) 


f iv-Vs)^ + iVs-Vi)^ \ 
\ Hvs - Vi)^ ) 


(77) 


This indicates that perturbations can grow before the singular point (coming from Regime II), and that they could 
grow after the singular point (toward Regime I), but we saw that this regime is not present in the bounce phase (Sec. 


VA), so the bounce phase will end shortly after the singular point rjg 


D. Discussion 


Let us summarize the evolution of curvature perturbations on super-Hubble scales through the bounce phase. 
Figure is a sketch of the evolution of C according to the results found above. If C enters the bounce phase with 
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a nonvanishing time derivative^, then we find that curvature perturbations can grow at most linearly in some time 
interval [tamp-, iamp+] and that the growth is maximal at the bounce point. This happens in what we call Regime II. 
We highlight this regime in purple in Fig. Regime III follows Regime II at which point blows up and becomes 
infinite at some time tg. At this point, curvature perturbations become constant, and the bounce phase ends shortly 
after. We highlight this regime in cyan in Fig. 

In the end, the amplification that C receives is dominated by the growth during the interval [tamp-)^amp+]- Between 
the beginning of the bounce phase and the beginning of the amplification phase, we expect little growth of the curvature 
perturbations, and so, the initial time derivative of Q at the beginning of the amplification phase should be of the same 
order as the time derivative of C at the beginning of the bounce phase, which we expect to be small. In fact, in the 
Ekpyrotic phase of contraction (where ru ^ 1) which precedes the bounce phase, the dominant mode of C, is constant 
in time while the second mode is decaying (as shown in Appendix]^. Hence, the amplitude of C at the end of the 
period of Ekpyrotic contraction is the same as the amplitude at the end of the matter phase of contraction (assuming 
for a moment that there is no intermediate radiation phase). Consequently, this could lead to a suppression of C(tB-)^ 
and hence to a suppression of the growth of C in the bounce phase since, as we argued, — C(^amp-)- 

The reason why we can match C and C at the end of the Ekpyrotic phase of contraction with the beginning of the 
bounce phase comes from the matching conditions of cosmological perturbations [sni^is^. These conditions impose 
that the gravitational potential $fc(ry) and the modified curvature perturbation variable Cfc(^) ni'e continuous across 
any transition (e.g. from the Ekpyrotic phase of contraction to the bounce phase). The variable Cfc is defined as [H] 

a^a+3C^(^) • (78) 


On super-Hubble scales {k ^ H), we note that the second term of the above expression is suppressed, so C.k — Cfc- 
Thus, Cfe must also be continuous across a transition. That is why the values of Ck and Cfc at the end of the Ekpyrotic 
phase of contraction are taken as the initial conditions of the bounce phase. 

At this point, we note that the maximal growth rate for C is given by 


Cmax — 



0(tamp— ) 


2 


(79) 


and that C grows at most linearly in time. Therefore, one can say that 

C(tamp+) ~ C(tamp-) ^ Ci^B-) ( “ ) (tamp+ ~ tamp-) ■ (80) 

\<P(tamp-) j 

Furthermore, since C receives essentially no amplification outside the interval [tamp-) tamp-i-]) we can place an upper 
bound on the total growth that curvature perturbations on super-Hubble scales receive from the bounce phase. 


^ C(tB-H) - C(tB-) ^ C(tB-) / 4'B \ . 

C(tB-) “ C(tB-) ^ C(tB-) ) 

where we divide the growth AC by the initial size of C before the bounce to get a dimensionless quantity. 


(81) 


E. Comparison with tensor modes 


We recall the equation of motion for 
reduces to 


tensor modes given by Eq. (12), which in the IR limit on super-Hubble scales 


h" + 2-h' = 0 . 
a 


(82) 


^ If enters the bounce phase with a vanishingly small time derivative, then curvature perturbations will remain constant throughout and 
exit the bounce phase unaffected. 
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Once again, we drop the k index when it is clear that the modes are in the IR limit. Close to the bounce point, 
we recall that the scale factor is almost constant, i.e. 0 ( 77 ) ~ as. Thus, we are left with the equation h" ~ 0, and 
consequently. 


Kv) - KVi) + h'{rii){'q - r]i ), (83) 

or, equivalently, 

h{t) h{ti) + - ti) . (84) 


Thus, as in the case of curvature fluctuations in Region II in the vicinity of the bounce point, there is a linearly 
growing mode. Dimensional analysis, however, tells us that this growing mode will not overwhelm the constant mode. 
The argument is as follows: we can estimate h{ti) to be of the order Mh{ti), where M is the mass scale at the bounce. 
On the other hand, we expect the time interval of the bounce phase to be of the order and hence we expect the 

linearly growing term to be comparable at the end of the bounce phase to the constant mode. 

Comparing the coefficients of the linearly growing modes of the curvature fluctuations and the tensor modes, i.e. Eq. 
(79) and the coefficient of the growing mode in Eq. (84), respectively, we see that it is the extra factor of [^B/0(<amp-)]^ 
in the coefficient of the scalar modes which leads to the enhancement of the scalar power spectrum relative to the 
tensor power spectrum. 


VI. A COMPREHENSIVE ANALYSIS OF THE PRODUCTION OF PRIMORDIAL 
NON-GAUSSIANITIES DURING THE BOUNCE PHASE 


Now that we have identified the conditions under which the tensor-to-scalar ratio can be suppressed, we turn to the 
study of how the bispectrum evolves during the bounce phase. We make use of the formalism developed in |49] (see 
also [671168]). 


Our starting point is the expression (36) for the three-point function. From this expression it is clear that the 


bispectrum builds up over time, which is to say that the three-point function after the bounce equals the three-point 
function before the bounce plus the result of integrating the right-hand side of (36) over the time interval of the 


bounce. From the form (38) of the interaction Lagrangian it follows that the terms which dominate the three-point 


function in the infrared are given by three powers of C and two powers of its time derivative. As shown explicitly in 
m in the computation of the three-point function in the matter-dominated contracting phase, the absolute amplitude 
of C cancels out in the definition of the shape function. Furthermore, Cai et al. |30| show that the bispectrum at the 
end of the period of matter contraction has an amplitude of the order 1 with a shape which is different from what 
is obtained in simple inflationary models. Since the dominant mode of ^ is constant during the Ekpyrotic phase of 
contraction, no additional contribution to the bispectrum is generated during that phase. We have not computed 
the contribution generated during a possible radiation phase of contraction between the end of the matter period 
and beginning of the Ekpyrotic period. This calculation could be done using the methods of |30j and we would find 
again a contribution with amplitude of the order of one and with a shape similar to that generated in the matter 
phase of contraction and different from that in simple inflationary models, the reason being that the same terms which 
dominate the bispectrum in the matter phase will also dominate in the radiation phase, and they are terms which are 
slow-roll suppressed during inflation. 

Hence, we now turn to the evaluation of the contribution of the bouncing phase to the three-point function. However, 
we must keep in mind that the equations of |49| , in particular the third order perturbed Lagrangian given by Eq. (38 1 , 


are only valid for a canonical scalar field. We must generalize the analysis to the case of the matter Lagrangian studied 
here (this generalization will not affect the evolution of the three-point function outside of the bounce phase because 
the extra terms which we derive below are negligible except in the bounce phase). This has already been done in the 
case of inflation for very general Lagrangians (see, e.g., [691 ED]). 

For the Lagrangian given by Eq. ([46]), perturbations up to third order in C yield the action 


53= J d'^x {Bi [dC.dxd'^C, - C,didj{d^C,djx)] + 


+ B^dx + + B.adCf + - 2/(C)^ ) , 


where 


^ ^ [acax - d-^d^dmd^x)] - 


(85) 


( 86 ) 
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The derivation of this action and the form of the functions An, Bn, and Cn {n = 1,...) can be found in Appendix 


Cl As expected, this action is equivalent to the action given by Eq. (38) in the limit where the Lagrangian (46) is 


canonical in a matter-dominated contracting universe. This is shown in Appendix |C 2 


In order to cancel the last term in Eq. (85), we make a field redefinition in Fourier space Q{rj, k) —>■ k) — f{r], k) 


in the third order Lagrangian. This way, there will be two contributions to the three-point function. The first part of 
the three-point function is the third order Lagrangian without the last term and the second part is related to the field 
redefinition terms where C(?7) k) is replaced by f{r], k). Using the Lagrangian formalism, we note that in Fourier space, 
we can canonically express the modes ({rj, k) as follows, 


C,{ri.k) = Cfc(??)4 + Cfc(??)a_ 


k ’ 


(87) 


where a^rjO) = 0, so is the annihilation operator, and at is the respective creation operator. Then, if we consider the 
interaction picture, the three-point function to leading order in the interaction coupling constant is given by 


{Civ, ki)C{v, k2)C{v, fc3))int =i df) ([C(J7, ki)C,{v, k2)C{v, ks), Lintiv)]) 


( 88 ) 


where rji corresponds to the initial time before which there is no non-Gaussianity. Also, Lint is associated with the 
third order action (85) without its last term. 

Here, we are interested in the production of non-Gaussianities during the bounce phase, so we consider the initial 
time to be the beginning of the bounce phase and we consider the end time at which the three-point function is 
evaluated to be the end of the bounce phase. However, as we saw in the previous section, curvature perturbations 
are nearly constant, and hence do not contribute to the three-point function, except during the small time interval 
[■'lamp-) ?7amp-i-] where C grows. Thus, the integration bounds are taken to be from ?7amp- to ryamp-i-j and the evolution 
of the curvature perturbations is taken to be 


Ckiv) = C'kivB-) + Ck'ivB-) 


4^' iV&mp—) 


iv V&mp-) ■ 


(89) 


The above expression follows from taking the maximal linear growth rate given by Eq. (79) throughout the amplification 


interval [r^amp-,?7amp+]- This expression slightly underestimates Cfc for ryamp- V Vb and slightly overestimates (k 
for rjB < V ^ f?amp-i- but it is a good approximation on average over the small interval [Tyamp-, ?7amp-i-]- 

We recall that curvature perturbations are more or less constant during the Ekpyrotic phase of contraction that 
precedes the bounce phase. Therefore, it is natural to take the end conditions of the matter-dominated phase of 
contraction as the initial conditions of the bounce phase. As shown in Sec. V D C,k and Cfc must be continuous across 


any transition on super-Hubble scales. Hence for the initial conditions of the bounce phase, we put the superscript 
which denotes the matter bounce solution EHI 


CTiv) = 


V2k^{r] - r)B-Y 


(90) 


where fjB- is the conformal time at the singularity if the matter-dominated contracting phase were to continue to 
arbitrary densities (i.e. without NEC violating matter). Also, A is a normalization constant which is determined from 
the quantum vacuum condition at Hubble radius crossing in the contracting phase, and it is found to be 


A = 


(AVB-r 

y/SasAfp 


(91) 


where At/b- = Vb- - Vb-- 

Let us comment on the wave number dependence of Eq. (89). We first solved the equation of motion in the bouncing 
phase in the limit where /c ^ 77 to 0*^ order. Then, matching the solution at the beginning of the bouncing phase with 
the one at the end of the matter contraction phase, we introduced some wave number dependence since the solution 
in the matter contraction phase has higher order terms in k/B. Thus, one may worry that obtaining the correct 
/c-dependent solution in the bounce phase up to leading order requires one to solve the full fc-dependent equation of 
motion. However, we note that we will be interested in the IR limit again when evaluating the three-point function. 
Thus, any k dependence not included in the above solution is suppressed during the bounce phase as long as k remains 
much smaller than the largest energy scale attained during the bounce, i.e. as long as k 'Hb-,Bb+, as long as 
the corresponding wavelength of the fluctuations remains much larger than the bounce length scale, which can be 
reformulated as fc ^ (Arys)”^. 
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Substituting the interaction Lagrangian Lint associated with the action (851 without its last term into Eq. ( 88 ) and 
using Eq. (87), we find 


(C(fci)C(4)C(4))i„t = (27r)35(3) CkAv+Kkiv+)Ck,{v+) 


rv+ 


xi dr] 


I h-h, 2 , [^2 • {h - hMh ■ {h + fcs)] 


M 2 


_ lu 

^3 




[ B2iv) j2 , -S3(77)^2(77) fci ■ fcs L 4 ( 7 ?)z^( 77 ) (fca ■ fca)^ -SrC??) 

i as ^ M'^aB k 1 M^as as 


+ 55 ( 77 ) 05(^1 • k 2 )CkAv)Ck 2 iv)Ck 3 {v) 


Ck^ivKkivKkiv) 


I CkAv)Ck2{v)Ck3iv) 

Ckiir])Ck2iv)Ck3iri) 

+ (5 permutations) . 

(92) 


Moreover, the contribution from the field redefinition is 

/ 3 \ 


-(C(fcl)C(fc2)/(fc3))redef = (2^)35(3) 

\i=l 

A2o(?7+)a| 


4M2 


2=1 / 

' t n t \ , (ki ■ k 3 )[{k 3 - iti) ■ ks]] 21 /- 7 m 2 

-fci • {ks - ki) + -p- I |Cfei(??+)l \Ck2{v+)\ 


Ais{r]+)aBz‘^{r]+) ( ki ■ (fca - fci) (fci • fc3)[(fc3 - fci) • fc3] 


M4 


kl 


^ 2^2 

^3 


/ 2A4(77+)a| - Ci(77+) 
V 2 z2(77+)c2 


Ck^{v+)Ck^(.V+)\Ck2{v+)\ 


CkSv+)CkM)\Ck2{v+W 

^ ' (5 permutations) . (93) 


The permutations that we refer to are over the ki vectors for i = 1,2,3. We note that, to simplify the notation, we set 
774 - = 77ainp± • The general form of the full three-point function can be expressed as 

(C(^l)C(fc2)C(4)) = (C(fcl)C(fc2)C(4))int + (C(^l)C(fc2)/(4))redef = (271)^5(3) ^ j fc 2 , fcs) , (94) 


and so, if we substitute Eqs. (|92[) and (93) into the above, we find the shape function to be given by 


A{ki,k2,k3) = 


klCkM)Ck2iv+)Ck2(v+) 

4Ci iv+)Q 2 (^+)Cfei iv+)Ck 2 iv+) 


rv+ 


X ^ 


dr] 


fri- 




kl ■ ks 2 [k2 ■ (fc 2 + fc3)][fc3 • (^2 + fcs)] 


_ 'hz _ 

fc 2 2 

fVQ 


fc2 


B2iv),2 , B3{ri)z‘^{r)) ki ■ k3 Bi{'q)z'^{'q) {k2 ■ k3f Bj{ri) 


Gb 


'fcf + 


M^ttB fc| 


M^qb 


7^2 i -2 
"-3 "-2 


Gb 


Ck^{v)Ck2{v)CLiv) 


Cfcl(?7)Cfc2(’7)Cfc3(??) 


55 ( 77 ) 05(4 • 4)Cfci(4Cfe2(4Cfe3(4 + ^^^^CkAA)Ck2iv)Ck3iv) 

®5 

^2o(??-H)a|fci / ^ (4 • 4 )[(4 - 4) ■ 4]^ 


4M2 


fc2 


24i8(77+)a5-z^(?7-i-)fc3 f 4 ■ (4- 4) (4 • 4)[(4-4) ■ 4]) Ckii'n+)Ckiiv+) 


M4 

iv±p 


fcl 


I 1,2 i-2t-2 

/2244(77+)a| -Ci{r]+)\ C'k,iv+)CkSv+) 


\CkAv+W 


\ 2z'^A+)cl 


KkAn+W 


(5 permutations) . 


(95) 


At this point, we should note that the contributions coming from the terms with coefficients Bi, B 2 , B 3 , and A 20 
are of order O(k^), and consequently, these terms are vanishingly small compared to other terms, which are of order 
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0{k^), on super-Hubble scales. Therefore, the three main contributions to the shape function are the term, the 
term, and the field redefinition term. We evaluate each of these terms separately in Appendix and we find the 
general expression for the shape function after the bounce phase [see Eq. (D18)]. 

Three important forms of non-Gaussianity in cosmological observations are the local form, the equilateral form, and 
the orthogonal form. The local form of non-Gaussianity requires that one of the three momentum modes exits the 
Hubble radius much earlier than the other two, i.e. ki k 2 = ks. Evaluating the shape function (D18l in this limit 
yields 


rlocal _ 10 A{ki ^2 — ks) 

/nl — o " 


10 

y 


-A 2 




B4{r]B)z% 


M4 


-b 


8aBAr]%_ V</''(?7-) 

Ais{'n+)aBz‘^{r]+) 2A4{r]+)a% - Ci{r]+) 


+ Bj{t]b) ) - 


3A2 


8 A 77 |_Ar 7 amp \4>'{ri^) 




Beivs) 


8M4A77ainp 8z2(7y_^)c2A?7amp 

The equilateral form of non-Gaussianity requires that ki = k 2 = ks, so one finds 


(96) 


^equil _ fO Ai^ki /C 2 ^ 3 ) 

/nl ~ ~Z, 


10 

y 


E.kf 

A^ 




B3{vb)z% ^4(?7b)4 


M 2 

-^p 


2M4 


- 257 ( 775 ) - 


3A2 


8Ar]%_Ari, 


amp 




BeivB) 


-b 


3Ais{r]+)aBz‘^{r)+) 2 A 4 ( 77 +)a| - < 71 ( 77 +) 


16M4A77amp 8 z2(77+)c 2A77 amp 

Finally, the orthogonal form of non-Gaussianity requires that fci = -y/-b = V^k, so one finds 


(97) 


^ortho „ 10 A{ki = y/kl + kl = V2k) 


/ortho _ 

/NL — 


Ey? 


10 

y 


A 2 


16aBA77|_ 
3A2 


8Ar]%_Ar]., 


amp 


(t^B \ ^ Bq(i]b) , q , f^^^ Ai^{A+)aBz‘^{il+) 2 ^ 4 ( 77 +)o| - (7 i(?7+) 

''{ri-)) a| 16M4A77anip 8z'^{r]+)c^^Ar]s,n,p 


(98) 


Substituting in some values for the model parameters (T, T, 7 , etc., introduced in Sec. IVI would yield specific 
numbers for the amount of non-Gaussianities that has been produced during the bounce. However, instead of giving 
exact values now, we will try to constrain the parameter space from observations. This is what we do in the next 
section. 


VII. COMBINATION OF THE OBSERVATIONAL BOUNDS ON NON-GAUSSIANITIES AND ON THE 

TENSOR-TO-SCALAR RATIO 


Let first rewrite the expression for using Eqs. (91), (D9), (DIO), (D12), (D16), and (D17l, 


/local 

/nl 


+ 37^'! + 12a|/3A/®T 


‘'B 
25apM^ 

W 

ba^B^P 1 


p 


37y'i 




Arj: 


amp 


•P'iv-] 


1 


A77amp 

IO 7 




1 


3//^5 A77amp 


16a|7"M4</>'jT) 

1 


47^05 


A77, 


(/)'(77_) 


</>'(7-) 


amp 


(/)'(77_) 


(99) 


The equilateral and orthogonal /nl have similar expressions, only with different coefficients. 
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At this point, we do not want to insert specific values for the model parameters. Yet, in order to have a healthy 
bounce, i.e. one that yields a bounce free of ghost instabilities, we expect the model parameters to he in specific 
regimes. From [UlllSlllH], we expect that (p'g < asMp, T /3 ~ 0(1), and 7 ^ 1. Also, from Eq. (54|, it is 

obvious that (j)'g/4>'[r]-) > 1. Therefore, keeping only dominant terms, the expression for reduces to 


/local 

/nl — 


5a%l3M^ f (j)‘ 




l(t>B 


3A77a 


+ 


1 


lAryamp Xfjj'ir]-) 


( 100 ) 


In the square bracket, the three terms come from Aqu, and Aredef, respectively. However, let us recall from 

Appendix [P] that the results for Acc^ and A(^ii were actually upper bounds in absolute value. Since we expect that 
Ar^amp (^l/asMp) from Eq. ( [70|, it results that the field redefinition term is dominant over the other ones, just 
like in the regular matter bounce |30|. so we can write 


/local ^ 

/nl — 


5ai/3M4 


A^amp V</>'(^-) 




( 101 ) 


In order to combine the bound on curvature perturbations and the above result, it is useful to rewrite the expressions 
for /nl in terms of the ratio A(/({t]b-)- In Sec. Eq. (81) told us that 


AC 


< 2 


CiVB-) ( <(’' 


C,{j]B-) ^ CiVB-) Wiv-) 


"B 


Arj; 


amp 


( 102 ) 


In the previous section, we argued that the initial conditions at rjB- were given by the end conditions of the 
matter-dominated phase of contraction, so we can say that 


(iVB-) CkiVB-) , 


CivB-)'^'^"\nVB-) ■ 


(103) 


Recalling that C™ is given by Eq. (90), we find that 

CiVB-) 


CiVB-) 


^ lim 

k/n^o Cr(^B-) 


3 

At]b- 


(104) 


and thus Eq. (102) becomes 


1 / ArjB- 

6 VA?7, 


amp 


\CiVB-) 


< 


(/'( 77 _) 


This allows us to place a lower bound on Eq. (101) as follows, 


/nl“‘ > 


54/^^^P 1 ( Ar^B- Y^ f AC 

144V6j3(/)'5 At] 

amp VA?7ampy \C{f]B-) 


5/2 


(105) 


(106) 


Our initial estimation in Sec. IIIC showed that we expected /nl to have terms of order (AC/C)^ and (AC/C)^- 
The terms of order (AC/C)^ in the full calculation corresponded to terms of order [(p'g/(p 'in our approximation 
scheme and they originated from Aqc,i 2 and A term of order (AC/C)^, i-e. of order [(f)'g/(f)', could have 


originated from Aqc,''^ but the full calculation showed that it did not have any real component [see Eq. (D5)]. Instead, 


the full calculation showed the presence of terms of order (AC/C)^^^, (AC/C)^^^, and (AC/C)^^^ coming from the field 
redefinition contribution to the shape function. In the large amplification limit, we are left with one term of order 
(AC/C)®^^ as shown in Eq. (106). 


Let us recall that A//C is bounded from below in order to satisfy the current observational bound on the tensor-to- 
scalar ratio r. Using the bound (31), we can further constrain the bound (106), 


/nl“‘ > 240 ( 


/ qbMp 


" ^ {ubMp)-^ 


AtJ: 


amp 


/ ArjB- \ 
\ A^amp / 


5/2 


(107) 


Let us note that Atjb- '^b- 

it results that ArjB- ^ {aBMp)~^. Thus, since every dimensionless ratio in Eq. (|107[) at least of order 1 or much 


'Hg_, and since the bounce energy scale is taken to be much less than the Planck scale. 


greater than 1, it results that /n“* > 240. Including the negative contribution to /n from the matter-dominated 
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contracting phase which is of order 1 |30| and the negative contributions from Ac^q'^ and Aqi3 would reduce this bound, 
but really not significantly. 

The best observational bounds on primordial non-Gaussianities currently come from the Planck experiment, which 
reports |32] 


/] 


local 

NL 


= 0.8 ±5.0 


/X‘' = -4±43 


f' 

J] 


ortho 

NL 


= -26 ±21 


(108) 


at 68 % CL. We see that the lower bound on coming from the bounce phase is already excluded by the observations 
at very high confidence level. Following the same analysis as above for the equilateral and orthogonal shapes yields the 
bounds > 359 and /nl^° ^ 289, which are also excluded at very high confidence level, although to a smaller 
extent than /n“*. 

To summarize, in this section we took the non-Gaussianity results derived in the previous section and imposed that 
there had been a sufficient amplification of curvature perturbations in order to satisfy the current observational bound 


on the tensor-to-scalar ratio. As a result, the theoretical lower bounds on /nl”* i /nl"° excluded at high 


confidence level by the current observational constraints on non-Gaussianities. This shows that the model suffers from 
the “no-go” theorem that we conjectured in Sec. HID 

Looking at Eq. ( flOTl ), we see that this could be alleviated if, for instance, the amplification period was very long 
compared to the Planck time, or if the model parameters were such that 1 or asMp/t/fg ^ 1. However, these 

conditions seem unlikely to occur in a physically admissible matter bounce scenario. 


VIII. CONCLUSIONS 

In the present paper, we have studied in detail the nonlinear dynamics of primordial curvature perturbations during 
the nonsingular bouncing phase in a matter bounce model described by a single generic scalar Held minimally coupled to 
Einstein gravity. This type of model can be made consistent with the observational bound on the tensor-to-scalar ratio 
by realizing an enhancement of the curvature perturbations in the bouncing phase. We derived the conditions on the 
model parameters for which such an enhancement can be achieved. We then expanded the action for perturbations up to 
the third order, computed the full set of three point correlation functions and then derived the nonlinearity parameters 
/nl in the cases of specific shapes of observational interest. Our results show that if the primordial non-Gaussianities 
mainly arise from a manifest growth of curvature perturbations during the bounce, then the nonlinearity parameter 
would become dangerously large and lead to a disagreement with the observational constraints from cosmic microwave 
background (CMB) data^. Specifically, we examined the relation between the nonlinearity parameter in the local, 
equilateral, and orthogonal limits and the growth of the curvature perturbations and explicitly showed that the 
observational bounds on the tensor-to-scalar ratio and the CMB bispectrum cannot be simultaneously satisfied. This 
leads us to conjecture that there is a “no-go” theorem for single field matter bounce cosmologies, assuming that the 
nonsingular bounce is realized by a generic scalar field minimally coupled to Einstein gravity, which would rule out a 
large class of matter bounce models. 

We note that this “no-go” theorem might be circumvented by dropping certain assumptions imposed above. For 
instance, if one introduces more than one degree of freedom such as in the matter bounce curvaton mechanism [SHI in], 
the constraints from the tensor-to-scalar ratio as well as from the primordial non-Gaussianities can be satisfied at 
the same time, the reason being that in the curvaton scenario the scalar fluctuations are generated by a different 
mechanism than the tensors. As another example, if the initial Bunch-Davies vacuum is noncanonical (e.g., in the 
ACDM bounce |6], the initial quantum vacuum has Cg ^ 1), the initial ratio of the tensor modes to the scalar modes 
can be suppressed, in which case there is no need for the curvature perturbations to be enhanced during the bounce. 

Our analysis also does not immediately apply to nonsingular matter bounce models in which the violation of the 
null energy condition is obtained by changes in the gravitational action (e.g., in Loop Quantum Cosmology izaizHi, 
Hofava-Lifshitz gravity [13], extended F{R) gravity [53l|54], modified Gauss-Bonnet gravity [55], or torsion gravity 
scenarios [MUST]). We might expect that the no-go theorem will extend to modified gravity matter bounce models 
which have the same number of degrees of freedom as General Relativity, in which case the tensor-to-scalar is generically 
large |28j . However, if the gravity model contains new degrees of freedom, then we might be in a situation similar 
to what happens in the curvaton scenario: the new degrees of freedom source the scalar modes but not the tensor 
modes, thus suppressing the tensor-to-scalar ratio during the bounce phase. Yet, it would be interesting to explicitly 


^ We recall that it has also been found in I.S.SI l,S4l that non-Gaussianities could become dangerously large in a certain nonsingular bouncing 
cosmology and it has been conjectured that this could be generic to a large family of nonsingular bouncing cosmologies. 
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analyze the conditions under which the bispectrum constraints can be made consistent with the observed bound on 
the tensor-to-scalar ratio in such models. 
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Appendix A: CURVATURE PERTURBATIONS EXPANDING ABOUT THE SINGULARITY 

The equation of motion for curvature perturbations in the IR limit is [see Eq. 0] 

dC , 


dr] 

Let us parametrize close to the singular point rjs as 


-c = 0 . 


so the equation of motion becomes 




dC 


iv - VsY 


<' ■ 


dt] rj-rjs 

Since after the singular time we have r] > rjs > ?7n we integrate as follows, 


In 


Civ) 

Civi) 


= 2 


dfj 


= 2 


V-Vs 


r^B — i r^s+f^ rri 


dfj 


+ej V-Vs 


(Al) 


(A2) 


(A3) 


(A4) 


for some constant e. As we take the limit e —>■ 0, the second integral vanishes by definition and we are left with the 
first and third integral. Evaluating them, we find 


In 


Civ) 

Civi) 


= 2 lim 
e-J-O 


In 


= 2 lim In 

e—i 

= 21n 


ivs -C-Vs 

Vi - Vs 

-eiv-Vs) 


-bln 


V-Vs 


iVs +C-Vs 


\ ( 77 * - r]s)e 
V-Vs 


Vs - Vi 


Therefore, 


Civ) = Civ^) 


V-Vs 


^Vs - Vi 

as expected if there were no singularity. A final integration yields 

' iv - Vs)^ + iVs - Vi)^ 


Civ) = Civ^) + Civ^) 


HVs - Vi)'^ 


(A5) 


(A6) 


(A7) 


As expected, we recover the constant mode solution (' = 0 as r] ^ 7]^. 
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Appendix B: PERTURBATIONS OUTSIDE THE BOUNCE PHASE 


Let us consider matter with an equation of state P = wp with w independent of time. In this case z{t) ^ a{t)Mp 
(see Sec. Ill A). Then, the solution to the long wavelength curvature perturbations is given by [see Eq. (72)] 


at) = au)+c{u)a{u)z{ur 
= au) + ma{uf 

Jti a (U 

For a constant w ^ — 1, the solution to the scale factor is given by 

a{t) = , 

for some positive constant oq, so we find 




au) + au)tr 
au) + ak)tf 



(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 


as long as jicj 7 ^ 1. Thus, for matter with jicj > 1, the solution for C exhibits a constant mode and a mode which is 
growing in an expanding universe (decaying in a contracting background), whereas for matter with jraj < 1 , it exhibits 
a constant mode and a mode which is decaying in an expanding universe (and growing in a contracting background). 
For example, this implies that an Ekpyrotic phase of contraction in which w ^ 1 has a constant mode and a decaying 
mode. 

For w = — 1, one would recover the standard result of inflation where the constant mode is dominant on super-Hubble 
scales in an expanding background, and the second mode dominates in a contracting space. 

The w = 1 case of fast roll expansion is relevant for the dynamics of our nonsingular bouncing cosmology right after 
the bounce phase. A phase of fast roll expansion occurs if the Lagrangian for the scalar field is dominated by its 
kinetic term, i.e. V{(j)) <^^/2. It then follows that the solution for the curvature perturbations in this case is (here in 

conformal time) 

({ 11 ) = C{V^) + [ -4^ ■ (B 6 ) 

Jr,, 

Solving the background dynamics tells us that the solution to the scale factor in a phase of fast roll expansion is 

a{v) = CEiv-VEy^"^ , (B7) 

where ce and pE are some constants. Thus, 

({v) = C{Vr) + CiVi){Vr-VE) [ (B8) 

Jr,, V-VE 

= C{Vr) + C{Vi){Vr-VE)^n(-—■ (B9) 

Km - veJ 

So, for w = 1, curvature perturbations exhibit a constant mode solution and a logarithmically growing mode, i.e. 
(■( 77 ) ~ In 77 . We note that this is also true in physical time since a ~ ~ 77 ^/^ implies that ^(t) ^ Int^/^ ^ Inf. 


Appendix C: THIRD ORDER PERTURBED ACTION 

1. Derivation of the general form of the third order action 

To study the three point correlation function in this matter bounce model, we have to evaluate the action up to 
third order in perturbation theory. We use the metric in the Arnowitt-Deser-Misner (ADM) form (see, e.g., [68]) 


ds^ = N'^dt^ — pijlN^dt + dx'‘){N^dt + dx^ 


(Cl) 
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where Ni = is the shift vector and N is the lapse function. The tensor 7 ^- is the metric of the 3-dimensional 

spacelike hypersurfaces in this 3 + 1 decomposition. It is related to the full 4-dimensional space-time metric tensor 
via \/—g = Ny^, where g and 7 are the determinants of the tensors g^^ and 7 ^-, respectively. The action (45) in this 
ADM decomposition is given by 


s= d^x 


^ - k" ) + X) + X)D(^ 


where is the three-dimensional Ricci scalar and the extrinsic curvature is defined by 

1 


— 2 jY DiNi DjNj) . 


We define the covariant derivative Di on the spacelike hypersurfaces such that it is torsion-free and satisfies 

= 0 . 

i(3) 

Then, R is the Riemann tensor associated with this connection, and 


n(3) _ 

^ ikj 


(3) 


i?(3) = 7 b R 

are the Ricci tensor and Ricci scalar, respectively. In the uniform field gauge where 

<50 = 0 , 

72 j — <r e ^ ; 

one can use the Hamiltonian and momentum constraints to determine the scalar contributions to the lapse function 
and shift vector. 


(C2) 

(C3) 

(C4) 

(C5) 

(C 6 ) 

(C7) 

(C 8 ) 


N = 1 + a , Nt = dia , 


(C9) 


up to leading order. Substituting Eq. (C9) into the metric [Eq. (Cl)] and expanding the action [Eq. (C2)] up to third 
order, we obtain the following. 


S 3 = / d'^x a3 


Q >2 >2 2 d(da 

aia + 02(0 + asCa + 04^-a + 05— —a + oeoCC + 


Mp {didjuY - {d'^crf 
2 


^dCdcr 


d'^c 




a + 2 M 2 rC^V - 2 ^ 200 ^ - - M^a 


OH 




(dcr 




.{dCy 


idCda 


- - 9M^C\ + 2MVAi^ + + 2M^CC—o -^ 


, <9^0- 


d^a Mp [didjaY — (i9^cr)^ 


-2M. 


2 di(djadidja 


p 


(CIO) 


where we defined the following coefficients, 


1 


11 


1 


ai = ^MpH — 0 ^ + q 4’ K,xxx ) — 2 iJ 0 ( 5G,x + —4> G,xx + -;4> G^xxx 

+ <^^ + —4'^G,x4> + q4>‘^G^^xx ) , 


a2 = - 9M2 r2 + 3^^ I iK^x + ^^"K^xx ) + 18i70" ( G,x + ^^'"G,xx ) - 30M G^^, + ^0"G,,^x 


1 


1 V 


1 > 


03 = — 6MpH + 6^^ ^G,x + ^ 0 ^G,xjf 

04 = 2MpH — 2<p^ ^G,x + ^0^G,xjf^ , 

05 = - 2M^H + 303g,x , 

06 = - 9 (-2M2i9 + 03 g,x) , 

07 = - 2M^H + 3<j)^G,x • 


1 
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We note that the Hamiltonian and momentum constraints yield (these can also be obtained by varying the action 
above with respect to a and a) 


respectively, where we defined 


2 A/ 2 C 

a = , 

u 

(Cll) 

5^(7 = ag(9\ + d'^x , 

(C12) 

u = 2M^H - Y^G^x , 

(C13) 

2 M 2 

Os = , 

(C14) 


and where 


^ “ M 2 ^ + n2 


+ 12ffrG,x + 3H^^G,xx - 2^0,^ - ^G^x^ 


(C15) 


If we substitute Eqs. (Clll and (C12) into the third order perturbed action [Eq. (CIO)], we obtain 

' 4 I 1 C" + ^2CC' + ^3(9'C)C' + A^idC)^ + A^dCdxC + AeCCd\ + AjCm" + A^dCOxC 


Sr = / d'^x 


a 

2 o 2 


+ H 9 C d x + AioC d C + ^iiC d X + Ai 2 diCdjCdidjC + Ai3diCdjCd,djX + Al^^iC^JX^^^JC 
+ Ai3di(^djxdidjX + (^leC + ^i7C)(^i^jC)^ + (^isC + Aig(){didjx)^ + ( 7 I 20 C + A2iC)didj(^didjX 


+ (7I22C + ^ 23 C)(^ C) + (^ 24 C + ^ 25 C)(^ x) + (^26C + ^ 27 C)^ C '9 X + ^28CC^ X 


+ Sf, 


(C16) 


where we defined the following, 
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^23 = TT^Os ! ^24 = - J i ^25 = ^ ^ Agp. = 


We note that Sb is a short-hand notation for all the boundary terms, which do not make a contribution to the 
calculation at 3rd order. After many integrations by part, we obtain 


Sr= J d'^x {Bi [dCdxd'^C, - C,didj{d^C,djx)] + H2C^5^C 


+ BsCdCdx + B^adrdjxY + B 5 C{ 5 C)" + i^eC" + i? 7 CX - 2 /(C) ^ ) , 
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where 
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(5/^2 2 AN 
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dC CL 


and where 
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2z^ci 


We also introduced the following, 
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Furthermore, 


C*! = A^a + (^27 — n42i) 


Mj 


and Cs is the speed of sound introduced earlier in Eq. (59). 


(C20) 


(C21) 


2. Third order perturbed action in the limit of the matter-dominated contracting phase 


Let us evaluate the third order action given by Eq. (C17) in a matter-dominated contracting phase when G((?i, X) = 0 
and K{(j),X) = M^X — V{(j)). In this case, we have 
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Thus, 




2a(t)H{tf 


B 2 {t) = 0 , 


Bi{t) = - 




= -2e{t)M^a{t) , 


H{tY 


2a{t)H(tY 2a{t) 


B^{t) = Mlait)- 
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Here, we consider = 2e(t)iJ(i)^, e(t) = — H(t)/H(t)'^, and Cg = 1. Therefore, we find that 
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This is equivalent to the third order action given in m noting that we defined d‘^x = whereas they considered 

d^X = C 


Appendix D: EVALUATING THE SHAPE FUNCTION IN THE BOUNCE PHASE 


We want to evaluate Eq. (95) which, as explained in the text, has three dominant terms: the term, the term, 
and the field redefinition term. Let us start with the contribution from the term to the shape function, which is 
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where we omit the 5 additional permutations for now. Using Eq. (89) for Ckiiv)^ 
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We note that taking . from Eq. (89) actually gives an upper bound on since Eq. (89) used the maximal growth 

rate (79) for the full range [r ]_, 77 +J. This introduces some small error in the final result but this will turn out to be 
unimportant since, as we will see, the field redefinition contribution to the shape function will dominate over this 
upper bound on . 4 <^c' 2 . 

The time-dependent terms that remain inside the integral are B 3 , B 4 , Bj, and (^ki- The latter, may 
experience a nontrivial enhancement during the interval [ 77 _,ry+], and consequently, it may contribute significantly to 
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the integral. The other terms, i.e. B^z^, and Sy defined in Eqs. (C20) and , contribute as almost constant 

terms in the integral over the range [r]-,ri+]. Therefore, we rewrite Eq. (D2| in the following form. 
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where we denote zb = z{rjB)- Performing the integral and using Eq. (891 for ^kiiv+) (again, this contributes to 
obtaining an upper bound for we obtain 
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where the modes are implicitly evaluated at tjb- and where we recall that 2A?7amp = V+ ~ V-- At this point, one 
could substitute C“(??b-) with Eq. (90) and write the full expression for Ac^c,'^. However, to satisfy the observational 
bound on the tensor-to-scalar ratio, we expect there to be a large amplification of curvature perturbations during the 
interval [7y_,77+]. In fact , from Eqs. (31) and (81), it must be that |C^'/CfcJI[<(’B/'/’^(^-)]^Ar?amp In that limit, 

the shape function (D4) reduces to (to leading order) 
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which is purely imaginary, and hence, does not physically contribute to the physical shape function. The next-to-leading 
order terms are 
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Using Eq. (90) for (y^{r]B-) and taking the limit k <^H, we find the leading order real-valued contribution to be 
H2 
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The BnirjB) terms can be evaluated using Eqs. (|C20[), ([7l|, and (55): 
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Similarly, we can use the previous procedure to find the contribution from the term to the shape function (again, 
omitting the additional perturbations for now), 
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The ellipsis in the third line denotes higher-order terms in |C™/C™'|[C '(^-) /CB]^(^^amp) and in the fourth line, we 
took the leading order real-valued term in the expansion. From Eq. ( |C20 ), the Bq term is given by 
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From the same argument as for the result (Dll) is actually an upper bound (in absolute value) for ^^/3, which 

we will comment on later. 

The contribution from the held redehnition term to the shape function is (again, omitting the additional perturbations 
for now) 
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where we took the leading order term in |^^/^^'|[(/)'(77_)/(/)^]^(A7yamp) Here, 
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Recalling Eq. (54) and the fact that [ 77 + — 77 b| = \r]- — tjbI = A? 7 amp, we have ((>'{'<1+) = 4>'{V-)^ so we can rewrite 
the terms above as 
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Combining the different contributions (including all permutations), the general form of the total shape function is 
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found to be 
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